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Results from unsteady Reynolds-averaged Navier-Stokes computations are described for two different syn- 
thetic jet flows issuing into a turbulent boundary layer crossflow through a circular orifice. In one case the jet 
effect is mostly contained within the boundary layer, while in the other case the jet effect extends beyond the 
boundary layer edge. Both cases have momentum flux ratios less than 2. Several numerical parameters are 
investigated, and some lessons learned regarding the CFD methods for computing these types of flow fields are 
outlined. Results in both cases are compared to experiment. 

I. Introduction 

Synthetic jets are attracting increasing interest in the aerospace community because they have many practical 
applications, including jet vectoring, control of separation, enhanced mixing, reduction of wall skin friction, and 
virtual aeroshaping. As a practical tool, they are generally more attractive than steady blowing or suction because 
they do not require any complex fluid ducting or plumbing systems. Furthermore, synthetic jets have been found to 
achieve similar effectiveness to steady blowing or suction with considerably smaller momentum (and hence lower 
energy inputs). 1 * 

The problem of synthetic jet flow in a turbulent boundary layer has been studied both experimentally (see, for 
example 2-5 ) and computationally (see, for example 6-11 ). Earlier numerical work by Rizzetta et al. 12 for synthetic jets 
into quiescent flow demonstrated that accounting for the internal actuator geometry was important because it affected 
the jet profiles at the exit. Lee and Goldstein 6 used direct numerical simulation (DNS) to study an array of synthetic 
microjets in a crossflow. A 4-slot array was modeled, and adjacent jets were coupled 180° out of phase. As microjets, 
the Reynolds number of the jets was very small, on the order of 10. Weak jets were confined to the linear sublayer, but 
when the jets penetrated up to the buffer region, they affected the events of the boundary layer significantly. 

Mittal et al. 7 and Mittal and Rampunggoon 8 used 2-D Navier-Stokes to investigate the use of stronger synthetic 
jets immersed in a boundary layer to form large mean recirculation bubbles (for aeroshaping). They modeled the 
cavity, but their studies were purely numerical (there was no associated experiment). Similarly, Ravi et al. 9 performed 
a purely numerical study of synthetic jet in crossflow. They also modeled the cavity, and used 3-D Navier-Stokes to 
investigate the effects of a square and two rectangular shaped slots. 

Cui et al. 10 performed 2-D computations on a synthetic jet in a turbulent crossflow boundary layer using un- 
steady Reynolds-averaged Navier-Stokes (URANS) with the Spalart-Allmaras turbulence model. 13 They compared 
their computational fluid dynamics (CFD) results with two experiments: one in which the discrete vortices scaled 
with the boundary layer thickness, and one in which they scaled with the inner viscous sublayer. Overall, they found 
qualitatively good agreement with experiment, although in the latter case the computed vortices moved too slowly 
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downstream. Also, in both cases there were significant quantitative differences between the CFD results and experi- 
ment. 

A CFD validation workshop on synthetic jets and turbulent separation control (CFDVAL2004) 14 ’ 15 was held in 
March 2004 in Williamsburg, VA. The purpose of the workshop was to assess the current capabilities of different 
classes of turbulent flow solution methodologies to predict flow fields induced by synthetic jets and separation control. 
Case 2 of the workshop was a synthetic jet issuing into a turbulent boundary layer through a circular orifice in the floor. 
This experiment was conducted at NASA Langley Research Center specifically for use as a CFD validation case. 16 A 
similar experiment was also recently conducted at NASA Glenn Research Center. 17 

A summary of the CFD results for case 2 of the CFDVAL2004 workshop is given in Rumsey et al. 15 There were 
five contributers who ran ten separate cases. One contributor used large eddy simulation (LES) and the others all used 
URANS. On the whole, reasonably good qualitative results were obtained compared with experiment both for long- 
time averaged and phase averaged quantities, but there were significant variations among the CFD results. Different 
turbulence models were found to have a much smaller effect than different grids, codes, and boundary conditions. Cui 
and Agarwal 11 also computed case 2 from CFDVAL2004, using both URANS with Menter’s SST turbulence model 18 
and detached eddy simulation (DES). 19 The two methods produced overall similar results in qualitative agreement 
with experiment, and DES was able to simulate some of the larger cross-stream velocity components downstream. 

The intent of the current paper is to summarize, compare, and contrast the CFD results obtained using the URANS 
code CFL3D 20 for the two cases hereafter described as the “CFDVAL case” (Schaeffler and Jenkins 16 ) and the “NASA 
Glenn case” (Milanovic et al. 17 ). Previous CFD results using CFL3D were reported for the CFDVAL workshop 
case, 21-22 but only limited results were reported for the NASA Glenn case in a paper devoted primarily to describing 
the experiment. 17 


II. Numerical Method 

The computer code CFL3D 20 solves the three-dimensional, time-dependent, Reynolds-averaged compressible 
Navier-Stokes equations with an upwind finite-volume formulation (it can also be exercised in two-dimensional mode 
of operation for 2-D cases). Because the code is compressible, the Navier-Stokes equations are averaged using Favre 
averaged variables; e.g., ~pui/~p. CFL3D can solve flows over multiple-zone grids that are connected in a one-to-one, 
patched, or overset manner, and can employ grid sequencing, multigrid, and local time stepping when accelerating 
convergence to steady state. Upwind-biased spatial differencing is used for the inviscid terms, and flux limiting is used 
to obtain smooth solutions in the vicinity of shock waves, when present. Viscous terms are centrally differenced, and 
cross-diffusion terms are neglected. For very low Mach number flows, preconditioning 23 is used to insure convergence 
and accuracy of the solutions. 

The CFL3D code is advanced in time with an implicit approximate factorization method. The implicit derivatives 
are written as spatially first-order accurate, which results in block tridiagonal inversions for each sweep. However, 
for solutions that utilize Roe flux-difference splitting, 24 the block tridiagonal inversions are further simplified using a 
diagonal algorithm with a spectral radius scaling of the viscous terms. 

In time-accurate mode, CFL3D uses pseudo-time stepping with multigrid and achieves second order temporal 
accuracy. With pseudo-time stepping, subiterations are used to reduce the linearization and factorization errors, and 
advance the solution in pseudo-time to the next physical time. For a non-deforming mesh, the time dependent com- 
pressible Navier-Stokes equations can be written as 


1 5Q 

J dt 


R{ Q) 


(i) 


where Q is the vector of conserved variables [p, pu, pv, pw, e] T , J is the Jacobian of the generalized coordinate 
transformation, and R is the right-hand-side: 


R{ Q) 


a(F-F„) d(G-G v ) d(H-H„) 

dt; dp dC, 


( 2 ) 


made up of gradients of the inviscid and viscous flux terms. Eq. (1) can be discretized with backward differencing: 


(1 + </)(Q” +1 - Q") - <fr( Q" - Q n_1 ) 
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( 3 ) 


The superscripts n,n + 1, and n — 1 indicate the time level. When <f> = 0 the method is first-order temporally accurate, 
and when 0 = 1/2 the method is second-order temporally accurate. To employ pseudo-time stepping, a pseudo time 
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term is added to Eq. (3), and the equation is iterated in m, where m is the subiteration counter. After linearizing R and 
performing some additional manipulation, one obtains: 
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where r is the pseudo time parameter, AQ'" = Q ,n+1 — Q m , and 


d(F - F v ) _ d(G - G v ) _ d(H - H„) 
dQ <9Q dQ 


(5) 


Eq. (4) is approximately factored and written in primitive variable form; it is solved as a series of sweeps in each 
coordinate direction. Additional details are given in Krist et al. 20 In the current study, the effects of different physical 
time steps and number of subiterations were investigated. 

The turbulence models are solved uncoupled from the mean flow equations using implicit approximate factoriza- 
tion. Their advective terms are solved using first-order upwind differencing. Many turbulence models are available 
in CFL3D, but only those used in the current study are mentioned here. Descriptions of the one-equation Spalart- 
Allmaras (SA) and the two-equation Menter k-ui shear-stress transport (SST) turbulence models can be found in their 
respective references. 13 - 18 Both of these models are linear eddy-viscosity models that make use of the Boussinesq 
eddy-viscosity hypothesis. The EASM-ko model is a nonlinear eddy-viscosity model: it stands for explicit algebraic 
stress model in k-u> form. A detailed description of this model can be found in Rumsey and Gatski. 25 

For time-accurate URANS computations, it can be shown that the dependent variables calculated for forced peri- 
odic cases correspond to phase-averaged variables. Following Reynolds and Hussain, 26 triple decomposition can be 
used to represent a fluctuating signal: 


/(x,t) = /(x) + /(x,f) + /'(x,f) (6) 

where / (x) is the (time) mean value, / (x, t) is the statistical contribution of the organized motion, and /' (x, t) is the 
turbulence. The long-time average is defined as 

/(x)= lim ^ [ f (x, t)dt (7) 

T—>oo 1 J Q 

and the phase average is defined as 


</(x,i)> = lim + nr) (8) 

N—>oo iV z — ' 
n—0 

where r in the case of a synthetic jet flow field is the period of the imposed oscillatory motion. Thus, the wave 
component is 

/( x ,£) = (/( x >i)) - /( x ) (9) 

and the decomposition of Eq. 6 can be written as 

/( x ,t) = (/(x,t)) + /' (x,t). (10) 

This form now resembles standard Reynolds splitting, except that the flow is split into a slowly varying mean 
(phase average) and a random fluctuating part. See also Gatski and Liu. 27 As a result, the final conservation equations 
in terms of phase-averaged variables ( ) are identical in form to the standard Reynolds-averaged equations. In other 
words, when URANS is used for a forced periodic flow field such as that produced by a synthetic jet, the computation 
typically eventually settles down to a nearly exactly repeatable periodic variation. If and when this repeatability is 
attained, any single point during the cycle corresponds to a phase-averaged result from the experiment. An example 
showing the cycle -to-cycle repeatability for a global integrated quantity is shown in Fig. 1 for the CFDVAL case. 
Here, coefficient of drag nondimensionalized by a unit reference length (the absolute level of which is not particularly 
meaningful in this case) is plotted as a function of iteration. An FFT analysis (not shown) over the last three periods 
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confirms that repeatability has been obtained, with the coefficients of lower frequency modes more than 2 orders of 
magnitude lower than that of the driving frequency. 

All of the results for the current two cases were expected to be symmetric about the center plane of the orifice. 
The CFDVAL case was computed with a full-plane grid, and results indeed maintained symmetry as expected. For the 
NASA Glenn case, use of a full-plane grid caused asymmetries to develop while establishing an initial steady-state 
flow field with constant blowing. The cause of these asymmetries is not precisely known, but it appears to be related 
to the strength of the blowing. Recall that the jet in the NASA Glenn case was strong enough to penetrate through 
the boundary layer to the freestream. If the jet strength was lowered significantly, then the full -plane constant blowing 
flow field remained symmetric. In any case, even with asymmetric initial conditions, these asymmetries went away 
and the solution became symmetric after the cyclic time-dependent boundary condition, to be described in the next 
section, was initiated. 


III. Results 

For the CFDVAL case, flow passed in and out of a circular orifice D = 6.35 mm in diameter. The orifice was 
located on the floor of a wind tunnel splitter plate with a turbulent boundary layer at M=0.1 (XJ 0 0 = 34.6 m/s) and 
approximate boundary layer thickness of 21 mm ( 5/D = 3.3). The Reynolds number was 2230 per mm, or 14,160 per 
orifice diameter. The jet was driven electro-mechanically by a bottom-mounted square-shaped rigid piston mounted 
on an elastic membrane inside the cavity chamber beneath the splitter plate. The cavity was approximately 1 .7 mm 
deep on average while the piston was operating, and the piston moved approximately ±0.77 mm. The thin volume 
inside the cavity varied during the operation of the piston, from roughly 1.0 x 10 4 - 2.6 x 10 4 mm 3 . The frequency 
was 150 Hz, and the maximum velocity (w max ) out of the orifice was approximately 1 .3//^. The momentum flux 
ratio, defined as (wmax/E^oc) 2 , was approximately 1.69. In this case, the effects due to the jet were mostly contained 
within the boundary layer. Experimental data were measured using laser-doppler velocimetry (LDV) and particle 
image velocimetry (PIV). 16 

For the NASA Glenn case, flow passed in and out of a circular orifice D = 19.05 mm in diameter. The Mach 
number of the freestream was M=0.0175 (U^ = 6.1 m/s), and the turbulent boundary layer thickness was approx- 
imately 7.6 mm {5/D = 0.4). The Reynolds number was 395.3 per mm, or 7530 per orifice diameter. The jet was 
driven by a loud speaker at the bottom of a cavity mounted beneath the tunnel floor. The cavity was 152 mm deep 
(with a volume of approximately 2.8 x 10 7 mm 3 ). The frequency was 24 Hz, the maximum velocity out of the orifice 
was approximately 1.0//^, and the momentum flux ratio was approximately 1.0. In this case the effects due to the jet 
extended well beyond the edge of the boundary layer. Experimental data were measured using hot-wire anemometry. 1 7 

For all of the results to follow, the .E-direction is downstream and ^-direction is up. The y = 0 plane represents 
the center plane aligned with the center of the orifice. Two different gridding approaches were taken for modeling the 
two cases. Because the CFDVAL case employed such a shallow cavity, the computational grid included the (as-built) 
cavity. A time-dependent velocity boundary condition was applied at the (fixed) bottom of the cavity to simulate the 
effect of a moving wall. Additional details can be found in Rumsey. 22 For the NASA Glenn case, the cavity was so 
large relative to the orifice that it was impractical to model it in the computations. Therefore, in this case the grid 
included the orifice only, which extended 25.4 mm below the tunnel floor. For both cases, the initial conditions were 
steady-state flow with constant blowing. Subsequently, a time-dependent velocity boundary condition was applied at 
the bottom of the orifice or cavity. This boundary condition was calibrated in order to attempt to match the velocity as 
a function of phase at the orifice exit, as measured experimentally. The time-dependent velocities prescribed were: 

u = 0 u = 0 w = [{pw) max / p\cos {2 t: Ft) (11) 

where ( pw) max = 0.0008/dooaoo for the CFDVAL case, and O.OlTSpooaoo for the NASA Glenn case. F was the 
frequency of oscillation. Due to the very low Mach number of the NASA Glenn case, preconditioning was employed 
in that particular computation. At the upstream grid boundary, a turbulent profile was specified that matched the 
corresponding experimental incoming boundary layer profile. 

Figure 2 shows time history of the velocities near the center of the orifice for the CFDVAL case. Results are shown 
for two different full-plane grids (the fine grid had 4.1 million total points and the medium grid was created using 
every other point from the fine grid), and for three different turbulence models (SA=Spalart-Allmaras, SST=Menter’s 
shear stress transport k-u>, and EASM=explicit algebraic stress model). For the most part differences due to grid and 
turbulence model were relatively small. The CFD mimicked the general overall trends measured in the experiment, but 
there were some noticeable specific differences. However, these results were the best it was possible to attain within 
the constraints imposed by the modeling and boundary conditions used. See Rumsey 22 for more discussion on these 
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differences. 

Figure 3 shows long-time average u- velocity contours for the CFDVAL case on the center plane using CFD with the 
SA model. Although not shown, results from the other turbulence models looked very similar. The mean perturbation 
to the velocity field was mostly contained within the boundary layer. There were not enough experimental data 
available to cover the same plotting range as Fig. 3, but Fig. 4 shows a comparison between CFD and experiment for 
(/-velocity contours in a zoomed-in region very near the orifice exit, where Xq is the orifice center. There was only a 
crude qualitative correspondence between the two. The CFD result predicted a larger mean perturbation to the flow 
field (including a very small reverse flow region near the orifice lip at (x — Xq)/D = 0.5 that is difficult to see in the 
figure), whereas the experimental PIV data showed the retarded velocity region to be smaller and offset downstream 
from the lip. Other results for the CFDVAL case can be found in Rumsey. 22 Generally CFD was successful in 
qualitatively predicting many flow features - both long-time-averaged and phase-averaged - compared with LDV 
measurements. 

For the NASA Glenn case, only one turbulence model, SA, was employed. As discussed in section II, there were 
some problems maintaining symmetry on a full-plane grid when establishing the constant blowing initial condition, 
but subsequent unsteady computations re-established symmetry. The effect of computing with a full-plane grid vs. a 
half-plane grid was investigated for this case. Results showed very little difference, so all NASA Glenn results for the 
remainder of the paper were obtained using the half-plane grid only. A sample comparison showing long-time average 
(/-velocity contours in a plane parallel to the floor is shown in Fig. 5. 

Figure 6 shows the fine half-plane grid, made up of 1.8 million points. This grid had the same (rescaled) distribution 
of grid points as the CFDVAL case outside of the orifice. A medium grid was also employed that was created using 
every other point from the fine grid. The time history of velocity at the orifice center is shown in Fig. 7. Two different 
grid sizes made very little difference in the results. The total velocity agreed fairly well with the experiment, which 
did not measure the individual velocity components like in the CFDVAL case. The largest discrepancy between CFD 
and experiment was in the magnitude during the suction part of the cycle. The computations accurately captured the 
double hump in total velocity during the discharge part of the cycle seen in the experiment. 

These results were obtained using a time step corresponding to 720 steps per cycle, and 10 subiterations per time 
step. The effects of time step and number of subiterations were also investigated for the NASA Glenn case. Total 
velocity at the orifice center is plotted in Fig. 8 for the three conditions of 720 steps per cycle and 10 subiterations, 
720 steps per cycle and 20 subiterations, and 360 steps per cycle and 10 subiterations. These results were obtained 
on the medium grid. There were imperceptible differences between the three results. Contour plots in the flow field 
also indicated a high degree of temporal convergence. For example. Fig. 9 shows phase-averaged //-velocity contours 
in the plane z = 0.192 mm at phase = 0°. The left half of the figure shows results using 360 steps per period and the 
right half shows results using 720 steps per period. 

Figure 10 shows long-time average contours of u/Uoo in the center plane for the NASA Glenn case. The agreement 
between CFD and experiment was qualitatively very good. The mean perturbation to the velocity field in this case 
extended well beyond the boundary layer edge. Phase averaged contour plots (not shown here, but given in Milanovic 
et al. 17 ) also showed good qualitative agreement. In particular, CFD captured the size and position of the disturbance 
very well as it convected down the plate after being emitted from the orifice, although the tilt in the structures predicted 
by the CFD was less pronounced than in the experiment. 

Quantitative comparisons of u- velocity profiles as a function of phase at a fixed position (x — xq) / D = 2.235, 
y = 0 are shown in Fig. 11. From Milanovic et al., 17 the convection speed of the synthetic jet structure created 
during the peak expulsion part of the cycle (phase ss 260°) was approximately 0.75(7 x . Thus, the structure convected 
to ( x — Xq )/D = 2.235 at approximately phase = 340°. As the figure shows, at the two phases 320° and 0° 
surrounding this “event,” the velocity profiles were significantly affected, and the CFD predicted perturbed velocities 
in the boundary layer that were too low compared with experiment. However, the profile shapes were qualitatively 
similar, and the profiles at the other three phases were quantitatively predicted very well. The long-time average 
also yielded excellent agreement. Comparisons of phase-averaged profiles for the CFDVAL case were reported in 
Rumsey. 22 Overall, the agreement for the CFDVAL case was generally better than for the NASA Glenn case during 
the “event” portion of the cycle. Figure 11 also compares results on both the fine and medium grids, giving an 
indication of the approximate effect of discretization error for these grid sizes. The effect was relatively small. 

Figures 12 and 13 are contour plots showing the vertical velocity at the orifice for both the CFDVAL and NASA 
Glenn cases, within 20° of peak discharge and peak suction. Note that for the two cases the phases (which were 
designated by the individual experiments) were not the same: for the CFDVAL case peak discharge occurred near 
phase = 120° and peak suction near phase = 300°, whereas for the NASA Glenn case peak discharge occurred near 
phase = 260° and peak suction near phase = 80°. However, this phase designation is arbitrary, and in the figure 
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the closest phase available from the stored CFD data is displayed. Both cases showed similar behavior. Near peak 
discharge, the velocities out of the orifice were very nonuniform, and the strongest velocities occurred near the side 
edges of the orifice, downstream of center. Near peak suction, the velocities were more uniform over much of the 
orifice area. These figures indicate the importance of modeling the inside of the orifice in CFD computations. Although 
a simple time-dependent top-hat profile might adequately model the relatively uniform inflow during the suction phase, 
it would be difficult to prescribe a time-dependent profile at the exit plane that mimicked the non-uniform velocities 
present during the expulsion part of the cycle. The effect of ignoring the inside of the orifice, and instead applying the 
time dependent boundary condition (Eq. 11) directly at the orifice exit, was investigated. Quantitative comparisons are 
given in Fig. 14, showing both phase-average results as well as the long-time average. Not modeling the inside of the 
orifice had an overall negative impact on the comparison with experiment, as expected. 


IV. Conclusions 

Time-dependent Reynolds-averaged Navier-Stokes computations were described for two different synthetic jet 
flows issuing into a turbulent boundary layer crossflow through a circular orifice. In one case the jet effect was mostly 
contained within the boundary layer, while in the other case the jet effect extended beyond the boundary layer edge. 
The former case was originally one of the subjects of the CFDVAL2004 workshop on synthetic jets. 

The use of URANS was shown to be appropriate and viable for synthetic jet flows. In both cases described here, 
the computations did a reasonably good job predicting the qualitative behavior of the synthetic jet flow field. However, 
quantitatively there were many specific differences. Several numerical parameters were investigated, including effect 
of grid size, time step, number of subiterations, and turbulence model. Differences between results on the fine and 
medium grids were relatively small. Using at least 360 time steps per period and at least 10 subiterations per time step 
was more than sufficient to produce results that did not change perceptibly with further temporal refinement. Different 
turbulence models for the CFDVAL case had some influence on the solution, but these differences were generally not 
too large. 

For synthetic jets like those considered here, either a full-plane grid or a half-plane grid (with symmetry) could be 
used. The current computations showed the importance of including the region inside the orifice when modeling the 
problem. Use of time-dependent top-hat boundary conditions at the orifice exit plane was an oversimplification that 
failed to capture the complex nature of the flow field near the orifice, particularly during the expulsion phase of the 
cycle. 
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Figure 1. Scaled drag coefficient as a function of iteration, showing periodic repeatability of URANS computation for CFDVAL case. 
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SA. fine grid 




Figure 2. Time history of velocity near center of orifice for CFDVAL case (z = 0.4 mm above the floor). 



Figure 3. Computed long-time-average u/U^ contours in center plane of CFDVAL case (orifice is centered at (x — xo) / D = 0), SA, fine 
grid. 
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Figure 4. Close-up of long-time-average u/Uoo contours in center plane of CFDVAL case near orifice; (a) CFD (SA, fine grid), (b) experi- 
ment (PIV). 



Figure 5. Contours of iong-time average //-velocity in the z = 0.192 mm plane for NASA Glenn case using SA, medium grids, showing 
difference between full-plane computations (left half) and half-plane with symmetry (right half). 
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Figure 6. Picture of the half-plane 6-zone grid used for NASA Glenn ease (1.8 million grid points total). 



Figure 7. Time history of velocity at center of orifice for NASA Glenn case, SA. 
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Figure 8. Effect of time step and number of subiterations on time history of total velocity at center of orifice for NASA Glenn case, SA, 
medium grid. 



left half: 360 steps/period right half 720 steps/period 


Figure 9. Contours of //-velocity at phase = 0° in the z = 0.192 mm plane for NASA Glenn case using SA, medium grid, showing difference 
between 360 steps per period (left half) and 720 steps per period (right half). 
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Figure 11. Profiles of u-velocity at (x — xo)/D = 2.235, y = 0 as a function of phase and long-time average for NASA Glenn case, SA. 


13 of 15 


American Institute of Aeronautics and Astronautics Paper 2005-5016 




(a) CFDVAL case (b) NASA Glenn case 

w/Uinf contours at phases 1 20, z-0 w/Uinf contours at phase-240, z=0 




Figure 12. Contours of ui-velocity at orifice exit plane near peak discharge part of cycle using SA, fine grid; CFDVAL case on left, NASA 
Glenn case on right. 


(a) CFDVAL case (b) NASA Glenn case 

w/UInf contours at phase=280, z=0 w/Uinf contours at phase=80, z=0 




Figure 13. Contours of ^/-velocity at orifice exit plane near peak suction part of cycle using SA, fine grid; CFDVAL case on left, NASA 
Glenn case on right. 
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(c) phase = 160 
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Figure 14. Effect of orifice treatment on phase-average and long-time-average u-velocity profiles at (a; — xo)/D = 2.235, y = 0 for NASA 
Glenn case, SA, medium grid. 
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